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q ; Abstract 

We derive a relation between four-fermion QED Green functions of different covariant gauges 
which defines the gauge dependence completely. We use the derived gauge dependence to 
check the gauge invariance of atom-like bound state calculations. We find that the existing 

t—i ' QED procedure does not provide gauge invariant binding energies. A way to a corrected 

^ ■ gauge invariant procedure is pointed out. 

<N ■ 

1 Introduction 

\o ■ 

QED gives a successful description of atom-like bound states. The recent measurement 
r£3 ' of the positronium life-time seems to remove the only discrepancy between theory and 
experiment in this field. Still, one can scrutinize general basis of the existing theory which 
involves far from trivial assumptions. The main one is that excited states correspond to 
simple poles of four-fermion QED Green function 0, [^, |J. In fact, one cannot prove it 
because of instability of excited states. Next, more technical, is that Bethe-Salpeter kernel 
is regular in total energy of fermions near the poles. Combination of the above assumptions 
leads to the generally accepted prescriptions (see, for example, @) for calculation of bound 



state parameters. Needless to say, any numerical success yielded by these rules supports but 
cannot prove the assumptions. 

Let us explain why it is doubtful that the above assumptions hold. To this end, consider 
propagator of a charged particle. Naively, one would expect that it has a simple pole at the 
particle mass. But it is well known (see, for example ||) that radiation of massless photons 
causes branch point singularity instead of the simple pole. One should expect the similar 
effect for atom-like bound states. The only difference is that two-particle bound state is a 
dipole. Consequently, one expects the radiation to be less important. This expectation is in 
accord with the successes of the standard approach to the atom-like bound states. 

The aim of this paper is to demonstrate that the main assumption — correspondence of 
excited states to simple poles of the Green function — is in contradiction with gauge invari- 
ance. More precisely, we will show that the assumption leads to a gauge dependence in the 
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pole positions, i.e., in observable energy level shifts. We will estimate the leading contribu- 
tion to the derivative of level shifts over gauge-fixing parameter. It will turn out that the 
gauge dependence is too weak to be seen in calculations performed up to now. 

It may seem that there is an opposite statement in the literature. Namely, it was found 
in that level shifts of the standard procedure are gauge invariant. The difference between 
H and the present paper lies in the assumptions on the Green function properties which 
were used in the study of gauge invariance. In fact, derivation of M is based on the above 
assumptions which we do not use in our analysis. To be specific, in the first of two papers 
U it was pointed out that derivatives of the Green function over gauge parameter contain 
branch points in the total energy of the pair. The authors conclude, seemingly using the 
assumption that the only relevant singularities are simple poles, that these branch points 
should be shifted from the poles corresponding to the bound states. In the present work 
we allow the possibility that the Green function have branch points and simple poles of 
coinciding positions. In the second paper of 0, an explicit form of the level shifts was used 
to prove their gauge invariance. The derivation is algebraic in nature and employs implicitly 
the second assumption — namely, that Bethe-Salpeter kernels and their energy derivatives are 
finite at the poles. (In notations of |J, that means finiteness of fcW, JfeW, jfeW, ...) A gam, we 
do not use any assumption on the Bethe-Salpeter kernels in our work (in fact, we even don't 
need these objects) and arrive at an opposite result. One may conclude that some of the 
quantities few, k^ l \ ... of f| are ill-defined. Indeed, more close analysis proves @, that, 
say, is infra-red divergent. We should stress that one would run into these singularities 
in the Bethe-Salpeter kernels only in a calculation of level shifts of order at . 

The latter may give a wrong impression that one can safely use the standard prescriptions 
for level shift calculations up to order a 11 . The real range of applicability of the standard 
prescriptions can be found only from a comparison with new, corrected prescriptions. We 
have not them in our possession. So, the only claim of the present paper is that the standard 
prescriptions break down in order a 11 . 

We should anticipate a question on the gauge dependence of the ground level shift which 
follows from our general formulas. Indeed, there is no doubt that ground level of bound sys- 
tem, if it exists, corresponds to a simple pole of the Green function. But since perturbations 
mix it with excited states, the lack of consistent picture for exited states causes inconsistency 
in its description as well. 

The last reservation we should make is on the dependence of the effect under consideration 
on the masses of bounded particles. To simplify the interim formulae, we consider only 
fermion-antifermion bound states. But all can be generalized for arbitrary mass ratio. The 
mass in the final formulae becomes then the reduced mass of the pair. Thus, we claim that 
even for the case of infinite mass of the heavier particle, i.e., when it can be replaced by 
the external Coulomb field, the effect survives. We expect that this case may be the most 
appropriate one to try to develop new, corrected prescriptions for level shift calculations. 

Turning to a description of the present work itself, its main technical means is an explicit 
form of gauge dependence of the four-fermion QED Green function. We found a relation be- 
tween the Green functions of different covariant gauges which defines the gauge dependence 
completely. The derivation is nonperturbative and the relation may present some interest 
in itself. It turns out that the gauge dependence has a simple form in the space-time repre- 
sentation. To use it, we formulate a procedure of extraction of level shifts form the Green 



2 



function in x-representation. Comparison of the gauge dependence of the Green function 
with the extraction procedure allows us to find the gauge dependence of the level shifts. We 
conclude pointing out a possible way to a corrected gauge invariant procedure of level shift 
calculations. 

Next section contains a derivation of the evolution in the gauge-fixing parameter; section 3 
comprises a brief recall of the extraction procedure and an utilization of the general evolution 
formula from section 2 for an analysis of gauge-dependence of the extraction; in the last, 
fourth, section we point out the reason for the gauge dependence and a way to the correct 
procedure. 



2 Evolution in Gauge-Fixing Parameter 

Let us consider the four-fermion QED Green function 

Gp(xf,Tf,Xi,Ti)=i / DipDA exp {iS QE D{ft)) {ip(x f)ip(x f))(ifj(xi)ip(xi)) 



where Xf(xj) is a coordinate of outgoing particle (antiparticle) and Xi(xi) is the same for 
ingoing pair. The definition of gauge- fixing parameter ft is given by corresponding photon 
propagator: 

D -w,*> = /c^(-*. + ^)i*-. (2) 

Our aim is to study the dependence of Gp on ft. To this end, it is useful to consider 
a Green function in external photon field, G(A), which is a result of integration over the 
fermion field in the rhs of eq. ([!]). From the one hand, it is simply connected to the Green 
function: 

Gp = (e L eG(A)) A=0 , Lp = l±-DMj^. (3) 

(In this formula each Lp generates a photon propagator; the dependence on the coordinates 
of ingoing and outgoing particles is suppressed for brevity.) From the other hand, G(A) is 
simply connected to a gauge invariant object G inv (A): 

G(A) = G mv (A) exp (ie J*' A^dx» - xe jf ' A^dx^j . (4) 

The gauge invariance of Gi nv means that it is independent of the longitudinal component of 
A: 

d^G inv (A) = (5) 
and is a consequence of gauge invariance of the combination 

ip(x)exp(ie J A^dz^j ip(y). (6) 
A substitution of eq. (f|) into eq. (Q) yields 

Gp = (e L 0G inv (A) exp lie H A ^ ~ ™ C A^dxA ) . (7) 

V V Jx f Jx i ) ) A=0 
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Let us take a /3-derivative of both sides of this equation: 



—G/s = \e L f ) {dpLp)G inv {A) exp |ie jT' A^dx^ - ie £ A^dx" 



A=0 



To get an evolution equation, one needs to express the rhs of this equation in terms of Gp. 
It is possible because (dpLp) commutes with G inv (A) and gives a c-factor when acts on the 
consequent exponential. So, eq. (g) transforms itself into 

where we have restored the x-dependence and used F to denote the action of (dpLp) on the 
exponential: 

(dpLp) exp (ie f_ f A^dx^ - ie (_ A^dx^ J = 

\ JXf JXi J 

F(xf,Xf,Xi,Xi) exp [ie / A^dx^ — ie (_ A^dx^j . (10) 

An explanation is in order: In deriving eq. (|Sp we have used a commutativity of (dpLp) 
and G im ,(74); it is a direct consequence of gauge invariance of G inv (see eq. @) and the 
fact that ippLp) contains only derivatives in longitudinal components of A (see eq. (|J) for a 
definition of Lp and eq. fl2|) for /3-dependence of -D M1/ ). 

The solution of eq. (§) for /^-evolution is 

Gp(x f ,Xf,Xi,Xi) = exp(((3 - Po)F(xf,Xf,Xi,Xi))Gp (x f ,Xf,Xi,Xi). (11) 



dpD^P, x) = -j^d^ ln((x 2 - ^)m 2 ), (12) 



To get the final answer one needs an explicit view of F from eq. (11) . It is easily deduced 
from the F-definition ( Jl0| ) and the following representation for the longitudinal part of the 
photon propagator: 

1 

16n 2 

where m is an arbitrary mass scale which is fixed, for definiteness, on the fermion mass. 
Then, up to an additive constant, 

F = ^( ln I | ln ( x f- x i?(xf-Xi) 2 \ (13) 

An \ m 4 (xf — x~f) 2 (xi — xA 2 (xj — Xi) 2 (xf — Xi) 2 J 

Substituting eq. (|i~3|) into eq. (|TID, we get our final answer for /3-evolution: 



Gp(x f ,Xf,Xi,Xi) 



Z{x f -x l f{x f -x i ) 2 TftP-AO 



X 



m 4 (x/ — Xf) 2 (xi — Xi) 2 (xf — Xi) 2 (xf — Xi) 2 
Gp (x f ,Xf,Xi,Xi). (14) 



The normalization Z is infinite before the ultraviolet renormalization. After the renormal- 
ization it is scheme-dependent and calculable order by order in perturbation theory. We will 
not need its value in what follows. 
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3 The Bound State Parameters And The Four-Fermion 
QED Green Function 

The four-fermion QED Green function contains too much information for one who just 
going to calculate bound-state parameters. One can throw away unnecessary information 
by putting center of mass space-time coordinate of ingoing pair and relative times of both 
ingoing and outgoing pairs to zero: 

G( et )p(t,x, r',r) = Gp(xf(t,x,r'),Xf(t,x,r'),Xi(r),Xi(r)) , (15) 

where the space-time coordinates depend on a space-time coordinate of the center of mass 
of the outgoing pair (t, x) and a relative space coordinate of outgoing (r') and ingoing (r) 
pair. In the case of equal masses 

r' _ r' 

Z/ = (t,X+~), Z/ = (t,X-~), 

^ = (0,|), 3?<=(0,~). (16) 

G(et)p still contains an unnecessary piece of information — the dependence on the center 
of mass space coordinate. The natural way to remove it is to go over to momentum represen- 
tation and put the center of mass momentum to zero. In coordinate representation, which 
is more convenient for gauge invariance check, we define the propagator Dp of the fermion 
pair: 

G (et)p (t, x, r', r) = D p (t, r', r)5(x) + . . . , (17) 

where dots denote terms with derivatives of <5(x). It is natural to consider Dp as a time 
dependent kernel of an operator acting on wave-functions of relative coordinate. In what 
follows we will not make difference between a kernel and the corresponding operator. The 
naturalness of the above definition of the propagator is apparent in the nonrelativistic ap- 
proximation: 

e i2mt Dp(t) « ]T#(t)e- i£o 'P(£ ), (18) 

Eo 

where the summation runs over the spectrum of nonrelativistic Coulomb problem and P(Eq) 
are the projectors onto corresponding subspaces of the nonrelativistic state space. One can 
obtain eq. (0) keeping leading term in a-expansion of the lhs if one will keep t oc I /a 2 and 
r',r oc 1/q (see [|, §]). The subscript on E is to denote that it will get radiative corrections 
(see below). The exponential in the lhs is to make a natural shift in energy zero. In what 
follows we will include the energy shift in the definition of Dp(t). 

The next step in calculation of radiative corrections to the energy levels is a crucial one: 
one should make an assumption about the general form of a deformation of the t-dependence 



of the rhs of eq. (18) caused by relativistic corrections. A natural guess and the one which 
leads to the generally accepted rules of calculation of the relativistic corrections to the energy 
eigenvalues (see, for example 0) is to suppose that one can contrive oscillating part of the 
exact propagator Dp from the rhs of eq. (|18|) just shifting energy levels and modifying the 
operator coefficients P(E ): 

Dp{t)= £ 6(t)e~< Eo+AE °) t Pp(E + A Eo ) + ..., (19) 
e +a Eq 
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where dots denote terms which are slowly-varying in time (the natural time-scale here is 
1/Eq). The additional subscript (3 on Pp is to denote that oscillating part of Dp(t) can 
acquire a gauge parameter dependence from relativistic corrections. 

The conjecture (|T9|) could be proven if the bound states were the eigenstates of the Hamil- 
tonian. But being unstable they are not. We will see that the conjecture (|T9| ) contradicts 
gauge invariance. Still it turns out extremely useful — the relativistic corrections calculated 
with it are in agreement with the experiment. Is it possible that another ansatz may be used 
instead of eq. fllTf ) preserving its advantage of success is an open question. 

Let us see how one can use eq. fllPP in energy level calculations. It is quite sufficient to 
consider Dp(t) on relatively short times when A^ t <C 1, E t ~ 1. For such times one can 
approximate Dp expanding the rhs of eq. (^) over A Eo t: 

Dp(t) ^J29(t)e- iE ^t k Af\E ), (20) 



where 



<(£o) = E ^^Pp{Eo + A Eo ). (21) 



An extraction of these objects from the perturbation theory is an interim step in the level shift 
calculations. (Here we should mention that in calculation practice A^\Eq) are extracted in 
momentum representation — i.e. not as coefficients near the powers of time but as the ones 
near the propagator-like singularities (E — E + ie)~( k+1 \) To come nearer to the level shift 
values, useful objects are 

f =Ef(^ ! - ( 22 ) 

Namely, as notations of eq. ( plf ) suggest, eigenvalues of should be equal to normaliza- 
tions of bound state wave functions which are driven from unit by relativistic corrections 
while the eigenvalues of Ap should be energy shifts to the k-th power times corresponding 
normalizations. Thus, the eigenvalues of 

^Ai k) +Ai k) u^r 1 

(23) 



AO) 

P 



b P ~ 2 
should be just energy shifts to the k-th power. Thus, we define 

S P = Sjf> (24) 

to be the energy shift operator: its eigenvalues are the energy level shifts caused by relativistic 
corrections. Our aim is now to check /3-independence of Sp eigenvalues. 

Some notes are in order: If the conjecture (|i~9"D is true Ap should commute with 
and the following relation should hold: 

Sp k) = [Sp] k . (25) 

We will use it in what follows. Another thing to note is that relativistic corrections affect 
the form of the scalar product of wave functions and, thus, one should add a definition of 
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operator products to the formal expressions (p3|),(25). But the level of accuracy to which we 
will operate permits us not to go into this complication and use the operator products as 
they are in the nonrelativistic approximation — i.e. as the convolution of the corresponding 
kernels. 

The way to the gauge invariance check of the energy shift calculations is clear now: 



Using the gauge evolution relation ( |14|) one should find the /5-dependence of Sp and then of 



its eigenvalues. As Sp is defined in eqs. (gD,d2|) through Af h s which are, in turn, defined 
in eq. (p0|) through the propagator Dp, the first step is to simplify eq. (|T4l) to the reduced 
case of zero relative time and total momentum of the fermion pair: 



Dp(t,r',r) 



V(4t 2 )) 



;i-((r' + r) 2 /(4t 2 )) 



x 



m 2 r' 2 m 2 r 2 



Dp (t,r',r). 



(26) 



The factor in the square brackets of the second line is time-independent and further factor- 
izable on factors depending on either ingoing or outgoing pair parameters. This reduce the 
influence of this factor to a change in the normalization of states. Being interested in gauge 
invariance of energy shifts, we omit this factor in what follows. Let us turn to the analysis 
of the influence of the factor in the first line of eq. (PB|). 

This factor is close to unit in the atomic scale r',r ~ 1/a, t ~ 1/a 2 . We will use its 
approximate form: 

(27) 



a 



r'r 



Factor « 1 + —(/?- /3 )— + 0(a 

27T t z 

One can read the dependence of Ag on f3 from eqs. P 



as 



A 



A (k) _ 



a (P~ Po) 
27i(k + l)(fc + 2 



rA (k+2) r 



(28) 



where r is the vector operator of relative position of interacting particles. The mixing of 
different ^4^'s with a change in the gauge parameter is due to the presence of 1/t 2 in the 



relations (EBT) and the fact that 



rhs of eq. (|27|). Finally, using the definition (|2 

A(°> w 1 

in the nonrelativistic approximation one can derive the following /3-dependence of Sp: 



(29) 



S:, 



s, 



(30) 



Treating the term in the last line of the rhs of the above relation as a perturbation, one can 
get an approximate value of the /3-dependent piece of the energy shift just averaging the 
perturbation with respect to the corresponding eigenstate of Sp . 

Thus, we get for the leading order of /5-derivative of an energy shift the following repre- 
sentation: 



d_ 



A, 



a 

'Yix V6 



\ (rSlr) - \ (S L rSlr) - \ «r^)) , (31) 
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where (...) means averaging with respect to the corresponding nonrelativistic eigenstate and 
the subscript L means the leading order in a-expansion. 

Eq. (|3T| ) is sufficient to define an order in a in which the energy shifts become gauge 
dependent: 

(W;^ (32) 

Here we have taken into account that r ~ 1/a and Sl ~ a 4 . 

To have a gauge dependence in any observable is clearly unacceptable. In the next section 
we will see how one should correct the above procedure of energy shift extraction from the 
QED Green function to get rid of the gauge dependence of energy shifts. 



4 A Way Out 

The procedure recalled in the previous section is based on the conjecture (|I9"D. A consequence 



of this conjecture is the gauge dependence of energy shifts of eq. ([31]) . One can conclude that 
the conjecture is wrong. In particular, as one can infer from eq. fl26f ), the operator coefficients 
near the oscillating exponentials in eq. ([19]) should get a time dependence from relativistic 
corrections. Even if in some gauge they are time independent, the gauge parameter evolution 
should generate a dependence which in the leading order in a reduce itself to the following 
replacement in eq. (fl9|) : 

P,(E + A Eo ) -> Pp(E + A Eo ) + (33) 

That S / 3(£'o) has nothing to do with energy shifts but will give contributions to Ac , (i£o)'s 
from eq. (^Cj). Being gauge dependent these contributions lead to the gauge dependence of 
energy shifts. 

The way to the correct procedure is to throw away terms like S ( g(£'o)/^ 2 prior to the 
definition of the energy shift operator. Thus, a necessary step in the process of extracting 
energy shifts from the QED Green function (and the one which necessity is not recognized 
in the standard procedure) is to calculate and subtract contributions like the last term in 
the rhs of eq. ([53]) from the propagator of the fermion pair. 

Below we report on a calculation of Hp(Eo) from eq. fl3~3"|). The most economical way 
to calculate it is to note that the energy dependence of the Fourier transform of the corre- 
sponding contribution to the propagator is 

(E-E )ln(-(E-E + ie)) (34) 

and that it comes from diagrams describing radiation and subsequent absorption of a soft 
photon with no change in the level Eq of the radiating and absorbing bound state. Similar 
contributions (with another power of energy before the log) are well known for the propagator 
of a charged fermion ||. 

It may be worth to note here that contribution of eq. fl34|) vanishes at E = Eq. This 
explains why such contributions are insignificant for practical calculations of the present day 
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accuracy. In particular, one can neglect them, despite the log-singularity, in the resonance 
scattering calculations and preserve the classic results of ||. 

The first step in our calculation is to present the pair propagator in the following form: 

D p (t) « (eS ierA MA w (i, ^) e_ierA(0) ) a=o ' (35) 

where L s is the same as in eq. (|3]) except a restriction on the momentum of photon propagator 
- the range of its variation is restricted to the soft region which border is of order of atomic 
binding energies; the exponentials with gauge potential are originated from the ones in eq. 
(0); D inv is a descendant of Gi nv from (0): to go over from Gi nv to D inv one should make all 
pairing of non-soft photons in Gi nv and all the reductions of space-time coordinates which 
was involved in going over from the Gp of eq. flip to the Dp of eq. (|T7|); at last, all gauge 
potentials in eq. (|35|) are taken at zero of space coordinate in accord with the <5(x) of eq. 
(|I7|). The difference between the lhs and the rhs of eq. (|35| ) does not contribute to the term 
under the calculation. 

The leading in the nonrelativistic approximation contribution to D inv is the same as 
for Dp — it is just the propagator of the nonrelativistic Coulomb problem. We explicitly 
calculate the leading contribution to the dependence of D inv (t,A) on the gauge potential 
in its expansion over soft momenta of the external photons. Not surprisingly, the dipole 
interaction of the pair with the external photon field arises in this approximation: 

D inv (t, A) « ^ - H c + er£(t)^J , (36) 

where H c is the Hamiltonian of the nonrelativistic Coulomb problem and £ is the strength 
of the electric field: 

£(t) = -A(t)+VA (t). (37) 



Substituting eq. ( |3Tf ) in eq. ( ff5| ) and keeping terms with only one soft photon propagator 
we get expressions which sum contains the term under calculation: 

e 2 (L s rA(t)D nr (t)rA(0)) A=0 , (38) 
e 2 (l s J dr 1 dr 2 D nr (t - T 1 )r£(r 1 )D nr ( y T 1 - r 2 )rf (r 2 )I>„r(r2)^ , (39) 

ie 2 (L s j ' dr(D nr (t-r)r£(r)D nr (r)T A(0) - (40) 

vA{t)D nr {t - r)v£{r)D nr {r))] , 

/ A=0 

where D nr (t) is the propagator of the nonrelativistic Coulomb problem from the rhs of eq. 
©• 

The next step is to pick out a contribution of a level E in eqs. (p8|),(p9|),(^0|). That is 
achievable by the replacement 

D nr {t) -> e- iEot 9{t)P(Eo). (41) 
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The last ingredient that one needs to calculate eqs. (p8D,(p9D,(^0|) is the time dependence 
of the soft photon propagators. It can be deduced from eq. (fj) as 



(L s Mti)Mt 2 )) = e (fa - t 2 f > tl) ^ _ ^ 

(L a A i (h)e j (t 2 )) = ^((ti-t 2 ) 2 >t c 2 



s 



(L&itjSjfo)) = e((t 1 -t 3 ) 2 >t 



27T 2 (t 1 -t 2 )3' 



'V 



(42) 



Here the ^-functions are to account for the softness of the participating photons (t c ~ 1/ £0 

.JO 



Taking eq. (|42|) into account we get the following contributions from eqs. 
(38) 
ED 



t 



71 



Le{t)e- iE ^-\p{E Q )vP{E Q )vP{E Q ) 

t 71 O 



}-e(t)e- lEot ®i (P(E )rP(E )r - rP(E )rP(E )) . 



71 



(43) 



The sum of the above terms yields the result of our calculation 

M E o) = 

t(P(E )rP(E )r - rP(E )rP(E )) ). 



-(lp(E )rP(E )rP(E ) + (-1 + ^)rP(£ )r + 
7r V3 2 



(44) 



One can explicitly check that /3-dependence of Hp{Eo) is the right one — i.e. if one 
subtracts the S-term from the propagator before the definition of the energy shift operator, 
the latter becomes gauge independent. Another observation is that the S-term cannot be 
killed by any choice of the gauge (in contrast to the case of charged fermion propagator 
where an analogous term is equal to zero in the Yennie gauge). 

Summing up, in this paper we derived a relation between QED Green functions of different 
gauges. We used it to check the gauge invariance of the energy shift operator. It turns out 
to be gauge dependent. This fact forced us to recognize that energy shifts are not one, and 
the only one, source for the positive powers of time near the oscillating exponentials in the 
propagator of the pair. We found a particular additional source of the positive powers of 
time which is responsible for the gauge dependence of the naive energy shift operator. We 
conclude with an observation that at the moment we have not a clear definition of the energy 
shift operator — to get it one needs a criterion for picking out contributions to the positive 
powers of time originating from the energy shifts. 

The author is grateful to A. Kataev, E. Kuraev, V. Kuzmin, A. Kuznetsov, S. Larin, 
Kh. Nirov, E. Remiddi, V. Rubakov, D. Son, P. Tinyakov for helpful discussions. This work 
was supported in part by Russian Foundation for Basic Research, project no. 94-02-14428. 
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